Control and manipulation of entanglement between 
two coupled qubits by fast pulses 



Ferdi Altintas and Resul Eryigit 

Department of Physics, Abant Izzet Baysal University, Bolu, 14280-Turkey. 
E-mail: altintas_f@ibu.edu.tr, resul@ibu.edu.tr 

Abstract. We have investigated the analytical and numerical dynamics of 
entanglement for two qubits that interact with each other via Heisenberg XXX-type 
interaction and subject to local time-specific external kick and Gaussian pulse-type 
magnetic fields in x-y plane. The qubits have been assumed to be initially prepared in 
different pure separable and maximally entangled states and the effect of the strength 
and the direction of external fast pulses on concurrence has been investigated. The 
carefully designed kick or pulse sequences are found to enable one to obtain constant 
long-lasting entanglement with desired magnitude. Moreover, the time ordering effects 
are found to be important in the creation and manipulation of entanglement by external 
fields. 
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1. Introduction 

Entanglement exhibited by composite correlated quantum systems has become one of 
the most widely investigated physics subject of recent years pQ. It provides promising 
resource for some quantum information tasks, such as quantum teleportation [2], 
cryptography [3] and computing jl]. For those applications, the generation, 
manipulation, detection and control of entanglement between the quantum systems 
should be precise and effectively done during carrying out of the required operations. A 
spin 1/2 particle in a magnetic field is found to be a suitable candidate as a qubit of a 
quantum computer [5], so the Heisenberg chains have been used to construct a quantum 
computer based on quantum dots [5], nuclear spins [6], electronic spins [7j and optical 
lattices jH]- Moreover, Wu et al. showed that one qubit gates in spin-based quantum 
computers can be constructed with a global magnetic field and controllable Heisenberg 
exchange interactions [9]. Imamoglu et. al. demonstrated that 1-D Heisenberg chain 
of spin 1/2 particles (qubits) can be used as a prototypical system to study the role of 
entanglement in quantum computational tasks |10j . 

Numerous studies have been devoted to the control, production and manipulation 
of entanglement in Heisenberg spin chains with the help of external magnetic fields [HI 
H21 [131 [HI UHl Uni HZl HH 111 ISOl [21]. Among them, Sadiek et al. investigated the 
time evolution of entanglement between two Heisenberg XYZ coupled qubits under an 
external sinusoidal time-dependent nonuniform magnetic field and they demonstrated 
that the entanglement of the system can be controlled and tuned by varying the 
time- dependent magnetic field and the Heisenberg exchange coupling parameters [12"]. 
Abliz et al. studied the entanglement dynamics of two-qubit Heisenberg XYZ model 
affected by population relaxation and subject to various types of external magnetic 
fields, such as an inhomogeneous static field, homogeneous exponentially varying and 
periodically varying magnetic fields, and they demonstrated that high entanglement 
can be produced, controlled and modulated with the help of external time-varying 
fields despite the existence of dissipation [2]. Huang and Kais demonstrated that 
the entanglement in one-dimensional spin system, modeled by the XY Hamiltonian 
can be localized between nearest-neighbor qubits for certain values of the step function 
external magnetic fields [TTj . In Ref. [TT], the possibility of creating and controlling 
entangled states by changing the relative phase of control pulses was investigated. 
Huang et al. demonstrated that the dynamics of the nearest-neighbor entanglement 
for one-dimensional spin system under sinusoidal magnetic fields displays a periodic 
structure with a period related to that of the magnetic field [IB]. Blaauboer and Di 
Vincenzo presented a scheme for detecting entanglement between electron spins in a 
double-quantum-dot nanostructure and they demonstrated separation and coherently 
rotating of entangled spins in quantum dots by using a time-dependent gate voltage 
and magnetic field [T9j . 

The form of the external magnetic field used in control of entanglement is an 
important parameter in manipulation of the entanglement between the qubits. As 
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summarized above, the most of the external fields studied for the control are in static, 
exponential, step or sinusoidal form. An external field of the kick or Gaussian pulse 
sequence form might provide a better control schema for entanglement. Along these 
lines Kaplan et al. and Shakov et al. have investigated the population and coherence 
dynamics of a single qubit under the influence of kick or Gaussian pulse sequence 
magnetic fields [221 E2J. It was shown that the instantaneous pulses (kicks) provide 
a full population transfer in a qubit from one state to the other which has a great deal 
of interest in quantum computing and control theory [23] . In our recent paper [20J, 
we have extended the schema introduced in Refs. [221 123] to two coupled qubits and 
have investigated the dynamics of entanglement between two Heisenberg qubits under 
the influence of strong delta function (kick) and Gaussian pulse type magnetic fields 
directed along the z-a.xes. The qubits have been assumed to be initially prepared in 
the separable 1 01) and the maximally entangled ^(|10) + 1 01 ) ) states. We have shown 
that the fast pulses provide an efficient way of controlling entanglement between the 
coupled qubits and also by this control it is possible to manipulate the transition from 
disentangled to the entangled states of the system. Longitudinal control fields were 
found to be ineffective in creating entanglement for the 1 1 1) and 1 00) type initial states 
or manipulating the entanglement of ^(|11) ± 1 00) ) type maximally entangled initial 
states. 

In the present study, we extend the work in Ref. [20] to consider the effect of the 
transverse field and analyze the dynamics of entanglement between two qubits that 
interact with each other via Heisenberg XXX type interaction that are subject to site- 
specific time-dependent magnetic control fields. Our aim is to widen the possibility of 
creating and manipulating entanglement with the tailored external kick and Gaussian 
pulse sequence type fields for arbitrary pure initial states. One important finding of the 
present work is the transverse external fields can control the dynamics of entanglement 
for two qubits for initial states that could not be manipulated with the longitudinal 
fields. We have also showed that the opposite scenario can take place for transverse 
external fast pulses. 

The paper is organized as follows. In Sec. |2j the model and basic formulation for 
the solution of time evolution are briefly discussed. In Sec. [3j Wootters concurrence as 
a measure of entanglement is introduced. Entanglement dynamics of two qubits under 
multiple kicks analytically and multiple Gaussian pulses numerically are investigated in 
Sees. [| and [5j respectively. In Sec. [6j we conclude with a summary of important results. 

2. The model and basic formulation 

In the present study, we consider two Heisenberg XXX-coupled qubits under time- 
dependent external magnetic fields in x-y plane. The time-dependent Hamiltonian for 
this model can be represented as [El |25] {Ti = 1): 

H(t) = H + H mt (t), (1) 
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where 

Hq = J Yl °\°2i 

i=x,y,z 

H in t{t) = i ^^(0(008(5)^ + sin(e)^). (2) 

A 2 = 1 

where a\ 2 (i = x,y,z) are the usual Pauli spin operators, Bi(t) and B 2 (t) are the 
external time-dependent magnetic fields acting on the qubits 1 and 2, respectively, J is 
the qubit-qubit interaction strength and 9 is the angle between the magnetic fields and 
the x-axes. For simplicity, we will assume < 9 < |. Here, we take J to be constant 
in time and assume all the time dependence in the systems' Hamiltonian H(t) comes 
only from H int {t). It is also possible to control qubit-qubit entanglement evolution by 
time dependent coupling strength instead of magnetic fields (26] . It was shown that such 
controls can be implemented physically, for example by using interacting flux qubits [26J . 

The most general form of an initial pure state of the two-qubit system may be 
given by the state vector |tf(0)) = ai(0) |11) + a 2 (0) |10) + a 3 (0) |01) + a 4 (0) |00), 

4 

where Oj(O) (i = 1,2,3,4) are complex numbers with ^|aj(0)| 2 = 1. Then, under 

Hamiltonian Q the probability amplitudes evolve in time according to Schrodinger 
equation as: 
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where B^(t) = ^e~ t0 B k (t) (k = 1,2) and Bkif)* is its complex conjugate. 
Mathematically and conceptually, it is convenient to write the formal solution of Eq. (|3]) 
in terms of the time evolution matrix U(t) as: 
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(4) 



Here all the time-dependence of the system is contained in the time evolution operator 
U(t), while the initial conditions are specified in |\&(0)). The time evolution operator 
U (t) may be expressed as: 

U(t) = t t~ l fo A ^ dt ' =T e~ { fo( 6 o+H int (t>))dt' 

= T E —h / H(t n )dt n ... / H(t 2 )dt 2 / H{t x )dt!. (5) 
^ n\ Jo Jo Jo 

Here T is the so-called the Dyson time ordering operator which arranges the operators 
in order of increasing of time [27], for example, T H{U)H{tj) = H{tj)H{U) if tj > U 
and H(ti)H(tj) otherwise. It gives rise to observable, nonlocal, time ordering effects if 
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and only if H(tj),H(ti) 7^ [22l E3]. Note that the non-trivial time dependence in 

•f— 

U(t) arises from time dependent H[t) and time ordering T- 
3. Measure of entanglement 

For two-qubit systems, as an entanglement measure, Wootters concurrence is a well- 
defined quantity (28] . Its value ranges from for a separable state to 1 for a maximally 
entangled (Bell) state. The concurrence function is defined as: 



C{p) = max{0, ^Ai - y/X 2 - ^A 3 - y/Xt}, (6) 

where A, (i = 1,2,3,4) are the eigenvalues of the matrix p(t)(a\ ® 62) p(t)* ® ^2) 
in descending order. Here p(t) is the density matrix of the qubits and p(t)* is its 
complex conjugate. According to Schrodinger equation ([3]), it should be noted that the 
probability amplitudes cij(t) (i = 1,2,3,4) evolve in time interdependently. From this 
point, it is convenient to consider a general time-dependent two-qubit pure state in the 
form: |tf(t)) = a x {t) |ll) + a 2 (t) |10) + a 3 (t) |01) +a 4 (t) |00). Then, it is straightforward 
to show that the concurrence function (JgJ) for the general pure state with density matrix 
p(t) = \*&(t)) (*&(t)\ is given by the simple equation as: 

C(p) = 2 \ ai (t)a A (t) - a 2 (t)a 3 (t)\ , (7) 

where the time-dependent coefficients may be given by Eq. Q as: 

a i {t) = y £U ij {t)a j {0), (8) 

where Uij(t) (i,j = 1,2,3,4) are the matrix elements of U(t). According to Eqs. (J?| 
and ([8j), to study entanglement dynamics of two qubits, one has to present the initial 
preparation of the qubits, i.e., eij(O), and the matrix elements of the time evolution 
matrix U(t). In this work, we will consider pure states as initial states. An initially 
pure state remains pure at all times under the dynamics given by Eq. Q. It is well 
known that for pure states entanglement can quantify all quantum correlations between 
two two-level systems, while such an identification is complicated for mixed states and 
entanglement signifies only a particular type of quantum correlation [29J. 



4. Entanglement dynamics of two coupled qubits under the influence of 
three positive kicks 

Depending on the complexity of the time evolution in Hi nt (t), the dynamics may 
or may not be solved analytically. Most of the studies employ numerical solutions 
of the time evolution in order to control entanglement between two level quantum 
systems [TTj, H21 [131 HH HS1 El HE1 EO]- However, analytic solutions are more convenient 
and easy to analyze, if they are available. For one qubit case, Shakov et al. listed some 
progressive approximations in which the time evolution may be solved analytically [23J . 
These approximations are the qubit having degenerate basis states, the external field 
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being constant or changed slowly, the field being oscillating with a frequency close 
to the resonance frequency of the energy splitting of a qubit so that RWA solutions 
exist, the field being perturbative and a sudden or series of fast pulses (kicks). The 
fast pulses are the most important limiting cases compared to the others because it 
was shown that they provide a full population transfer from one state to the other in a 
qubit [221 [23] . Moreover, the fast pulse and pulses have many potential implementations 
for some quantum tasks, such as NMR, quantum gates, excitation of electronic states 
in molecules, chemical reactions and quantum computing [301 ED E21 E31 [311 EH] • 

We have obtained analytical expressions for the time evolution operator for one, 
two and three kick sequences and display the elements of U (t) in Appendices A, B 
and C, respectively They will be used in Eqs. and ^ to calculate the dynamics of 
concurrence for various initial states. In Ref. [20], it was found that the entanglement 
of the qubits initially in |01) (or |10)) and ^=([10) ± |01>) can be manipulated easily 
with an external field in ^-direction, while for the initial states of the type |11), 1 00) or 
linear combination of them were immune to such a manipulation. So, we consider those 
states that cannot be affected by longitudinal external fields and analyze their dynamics 
under a transverse field. 

In the following, we will discuss the effect of transverse field on the dynamics of the 
system which is initially in a separable, i.e., 1^(0)) = |11) or in a maximally entangled 
state i.e., |tf(0)) = ^(|H) + |00». To see the effect of kicks on the entanglement 
dynamics for these initial states, one should note that before the kick the qubits evolve 
in accordance with the time independent Hamiltonian Hq and the propagator is given 
by 

-iH t 

e~ iJt 

e iJt cos(2Jt) -ie iJt sin(2Jt) 
-ie m sin(2 Jt) e iJt cos(2Jt) 
e~ iJt 

Based on Eqs. ([7]), ([8f and ffijb, the initially separable state remains separable, while the 
concurrence for the initial Bell state is equal to 1 at any time before the kick. That is, 
qubit-qubit interaction has no effect in the absence of external fields. 

Here we consider two qubits whose states are strongly perturbed by three positive 
kicks applied at times t = T\, t = T 2 and t = T 3 . The magnetic fields on qubits can be 

3 3 

given as Bi(t) = a^5(i — 2$) and B 2 (t) = /3 — Tj), where a and (3 are called 

i=l i=l 

integrated magnetic strengths [20J, and 5(t) is the dirac delta function. For such a kick 
sequence, the integration over the time is trivial and the time evolution matrix ^ in 
the presence of time ordering can be obtained easily [201 [221 123]. I n the Appendix part, 
the evolution matrices after each kick are presented and by using these propagators 
in Figs. 1(a) and 1(b), the effects of the sequence of three kicks on C(p) for |11) and 
^(|11) + 1 00)) initial states are displayed, respectively, with 6 = n/2 and a/ (3 = 2 



U(t) 



(9) 
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Figure 1. C(p) versus Jt for a sequence of three positive kicks applied at times 
Ti = 5,T 2 = 10 and T 3 = 15 for the initially pure states |11) (a) and 4=(|11) + |00)) (b) 
with 9 = ^ . Here the dashed lines correspond to a/ ft = 2 and the solid lines to a/ ft = 3, 
and we set ft = 1 and J = 1. 



or a//3 = 3. The concurrence for the considered initial states in the time domains 
5 < t < 10, 10 < t < 15 and t > 15 may be found by using Eqs. ( ]A.2| ), ( |B.2[ ) and ( |C.2[ ), 
respectively, in Eqs. ([7]) and (|8]). Both figures show the pronounced effect of the kick on 
entanglement dynamics; the concurrence of these initial states starts oscillating just after 
a single kick with an increase in their oscillation amplitudes with the increase in the ratio 
a/(3 (see the time domain 5 < t < 10). In fact, a sudden kick applied at t — 5 induces 
entanglement from the initially unentangled qubits (Fig. 1(a)) and yields oscillations 
in the concurrence function for the Bell state (Fig. 1(b)). It is worth mentioning here 
that for the magnetic fields acting in the z-axes, it was found that the qubits remain 
unentangled or maximally entangled state at all times for |11) or ^(|11) + 1 00) ) initial 
states despite the existence of strong external kicks [20]. Comparing one, two and three 
kick regions, the sole effect of the number of kicks is found to change the amplitude (or 
minimum) of C(p) for initial Bell state. As can be seen from Fig. 1(b), after second 
kick, the amplitude of concurrence oscillations increases compared to that of first kick, 
while it decreases after third kick. On the other hand, for |11) state case, the effect 
of the number of kicks is complicated compared to the Bell state case, because each 
kick can change the amplitude, maximum and minimum of C(p). It is obvious that 
after each kick the minimum of C(p) increases for initially separable state. The most 
important observation from these figures is the possibility of obtaining almost steady 
high concurrence around 1 after third kick for a//3 = 3 and both initial states. This 
shows that for certain system parameters (here we set a — 3, /3 = 1, J — 1 and 
T = 5), it is possible to create highly entangled qubits from an initially separable state 
by perturbing the qubits via instantaneous pulses. Also note that after second kick, 
the minimum of C{p) for Bell state can go to for a/fl = 2. As mentioned above, 
the qubit-qubit interaction has no effect on the entanglement of the considered states 
before the kick, while their oscillation frequencies depend on the qubit-qubit interaction 
strength, J, after the kick as can be seen from Figs. 1(a) and 1(b) as well as from the 
matrix elements in the Appendix part. 




Figure 2. (Colour online) (a) shows the contour plot of C(p) versus Jt and a/ (3 with 
9 = ^. (b) shows the contour plot of C(p) versus Jt and 8/tt with a/j3 = 3. Here 
the contour plots are for |11) initial state and include three positive kicks applied at 
times T\ = 5,T2 = 10 and T3 = 15, and we set j3 — 1 and J = 1. In these contour 
plots we have assumed twenty equidistant contours of concurrence between (black) 
and 1 (white). 

To further elucidate the effect of the strength and the direction of the magnetic 
fields on C(p), we have displayed the contour plot of C(p) versus Jt and a/ {3 for 
the separable initial state in Fig. 2(a) and Bell state in Fig. 3(a) with 8 = 7t/2, 
while in Figs. 2(b) and 3(b), C(p) versus Jt and 8/tt have been plotted for |11) and 
^(|11) + 1 00) ) with a/ (3 = 3, respectively. In these figures, twenty equidistant contours 
of concurrence between (black) and 1 (white) are shown. From Fig. 2(b), it can be 
deduced that the entanglement dynamics for the initially separable state is independent 
of the direction of external fields as long as the field is transverse, which also can be seen 
from the relevant time evolution matrix elements in Appendices where the concurrence 
is ^-independent. Fig. 2(a) shows two pronounced results for a//3-dependence of the 
concurrence for |11) state: i) The separable state is found to remain almost separable 
despite the strong external kicks for a/ (3 = 1 and a/ (3 = 7.25, 13.25, 19.5. ii) The long 
lived high entanglement regions which are indicated in white straight stripes sections 
have long lifetimes only after the second and the third kicks. From the comparison of 
the one, two and three positive kick sections, it seems that each kick widens the long 
lived high entanglement a//3-area. This shows the necessity of using multiple kicks to 
create almost steady high entanglement. Also note that each kick creates a different 
a/ (3 oscillatory structure for C(p). 

Contrary to the initially separable state case, the concurrence for the initial Bell 
state strongly depends on the direction of the external magnetic fields if 8/tt > 0.04 
as can be seen from Fig. 3(b); for 8/tt < 0.04, the entanglement between two qubits 
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Figure 3. (Colour online) (a) shows the contour plot of C(p) versus Jt and a/ (3 with 
9 = ^. (b) shows the contour plot of C(p) versus Jt and 8/ir with a/ (3 = 3. Here 
the contour plots are for ^(|11) + 1 00) ) initial state and include three positive kicks 
applied at times Ti = 5,T 2 = 10 and T 3 = 15, and we set /3 = 1 and J = 1. In these 
contour plots we have assumed twenty equidistant contours of concurrence between 
(black) and 1 (white). 



is almost unperturbed by the external positive kicks. Similar to |11) state case, the 
concurrence for ^(|11) + 1 00) ) is found to be unaffected from three kicks for a//3 = 1 
and a/(3 ^ 4.25, 7.25, 10.5, 13.5, 16.75, 19.75 that can be deduced from Fig. 3(a). This 
figure also shows that every kick region has a different a//3 periodic structure and 
long-lived high entanglement can be effectively obtained and controlled only after third 
positive kick for the initial Bell state after perturbing the entanglement dynamics 
of qubits with external kicks, for example in the time domain 15 < t < 20 for 
a/0 = 2.75,5.25,9.25,11.75,15.5,18.0. From Fig. 3(b), one can see that it is possible 
to get a constant long-lived entanglement with desired magnitude after third kick by 
adjusting 9. For example, C(p) « 0.93, C{p) « 0.44, C(p) « 0.08 and C(p) « 0.63 for 
0/ti = 0.06, 6 /it = 0.18, 6/tt = 0.25 and 6/tt = 0.36, respectively. 

5. Entanglement dynamics of two qubits under the influence of Gaussian 
pulses 

The kicked approximation is based on the energy level of the qubit, AE, and the width 
of the pulse, r, and is valid if AEt << 1. Depending on the physical implementation of 
the qubit, it might be difficult to obtain an external field as a delta function kick. Instead 
of kicks, Gaussian pulses can be used. As is well known, the kicked approximation is 

/ (t-T K ) 2 \ 

the limiting case of Gaussian pulses I i.e., lim^o 3§^ e T ' 2 = &Ko{t — T K ) J , thus 
Gaussian pulses enable us to consider the effects of using finite-width pulses on 
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entanglement dynamics. From this point, in this section, we present the results of 
numerical calculations of concurrence for two qubits whose states are strongly perturbed 
by series of three positive narrow Gaussian pulses for the initially separable |11) and 
maximally entangled ^(|11) + 1 00) ) states. The concurrence for the considered initial 
states may be calculated by solving numerically the set of first order coupled ordinary 
differential equations in Eq. (J3|: 

ia^t) = J ai (t) + ~e- i6 B 2 (t)a 2 (t) + ~e~ ie B x {t)az{t) , 

ia 2 (t) = -e i6 B 2 (t)ax{t) - Ja 2 (t) + 2Ja 3 (t) + '£i(t)a 4 (t), 

ia 3 (t) = ~e*B 1 {t)a 1 (t) + 2Ja 2 {t) - Ja 3 (t) + ^e' 18 B 2 (t)a 4 {t), 

za 4 (*) = ^e te B 1 (t)a 2 (t) + l -e i6 B 2 (t)a 3 (t) + Ja 4 (t), (10) 

Jt, _(t-Tj) 2 - Jt, (t-Tj) 2 

with replacing B x (t) -»■ ^7=7 X, e r2 anc ^ B 2 {t) ->■ ^=v2_, e ^ . Here the 

i=l i=l 

Gaussian pulses are assumed to have the same width r and centered at times t = T±, T 2 
and T3. We will determine how the concurrences of the initially separable and Bell 
states depend on the pulse width. It should be noted here that if r were chosen small 
enough, the results obtained in the previous section could be reached. 

In Figs. 4(a)-4(f), we have shown C(p) versus Jt for a system strongly perturbed 
by three narrow Gaussian pulses for |11) and ^(|11) + 1 00) ) initial states with 9 = | 
and a/(3 = 2 or a//3 = 3. Comparing Figs. 1(a) and 1(b) with the subfigures in Fig. 4, 
the width of the pulse changes the minimum, maximum and the amplitude of C (p) for 
1 11) state, while for Bell state, this effect is in its minimum (or amplitude). The most 
pronounced observation from these figures is the existence of constant high concurrence 
nearly 1 at times t > 15 for a//3 = 3 and Jt = 0.3 for both initial states as can be seen 
from the solid lines in Figs. 4(e) and 4(f). From Fig. 4(e), it is safe to deduce that with 
a sufficiently high pulse width and a//3 ratio, it is possible to obtain a state very close 
to Bell state from an initially separable state by using a Gaussian pulse sequence. 

To further elucidate the interplay between a pulse sequence, relative magnetic 
strength on qubits, pulse width and entanglement, we display the contour plot of C(p) 
as function of Jt and a//3 at three different pulse widths Jt = 0.1, 0.2, 0.3 in Figs. 5(a)- 
5(c). Comparing the subfigures in Fig. 5 with the ideal kick case (Fig. 2(a)), the most 
important effect of the increasing pulse width is found to be enlargement of the area 
of the long lived high entanglement regions which are the white straight areas toward 
the top of Figs. 5(a)-5(c). For the kick case, it was concluded that the long lived 
high entanglement regions only appear after second and third positive kicks, while for 
wider Gaussian pulses, these regions also appear after the first pulse (for example, 
see the region 5 < t < 10 in Fig. 5(c)). Contrary to Fig. 2(a), a/ (3 periodicity does 
not exist in Figs. 5(a)-5(c); increasing the pulse width coalesces the a//3-dependent 
oscillatory structure and produce non-periodic structures. Furthermore, |11) state 
remains separable if and only if a/j3 — 1 under the strong influence of Gaussian pulses. 
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Figure 4. C(p) versus Jt for a sequence of three positive Gaussian pulses centered 
at T\ — 5,T 2 = 10 and T 3 = 15 for the initially pure states |11) ((a), (c) and (e)) and 
^(|H) + 1 00)) ((b), (d) and (f)) with 9 = § . Here the dashed lines correspond to 
a//3 — 2 and the solid lines to a//3 = 3, and we consider three dimensionless pulse 
widths as Jt = 0.1 ((a) and (b)), Jt = 0.2 ((c) and (d)) and Jt = 0.3 ((c) and (f)), 
and we set J = 1 and /? = 1. 



On the other hand, when r — > (see Fig. 2(a)), this holds for a//3 = 1 as well as for 
a/ (3 = 7.25, 13.25, 19.5, as was the case for the kick sequence. Peculiarly, Fig. 5(c) shows 
that there is a sudden transition between long lasting maximally entangled and almost 
separable states of the system just after a pulse for a = 13(3 and a = 19.8/3. Moreover, 
the subfigures in Fig. 5 demonstrates the possibility of obtaining desired values of steady 
entanglement by carefully designed pulse or pulse sequence and system parameters. 

In Fig. 6, the contour plot of C(p) versus Jt and a/ (3 are displayed for ^(|11) + |00)) 
initial state for the system under three Gaussian pulses having dimensionless pulse 
widths as Jt = 0.1,0.2 and 0.3. Similar to |11) state case, C(p) of Bell state is 
undisturbed by the highly wider Gaussian pulses if and only if a/ (3 = 1. In Fig. 3(a), 
we have noted that the long lived high entanglement regions can only be obtained after 
third kick after disturbing the entanglement dynamics for the Bell state, while these 
regions also appears after applying first and second Gaussian pulses as can be seen from 
Figs. 6(b) and 6(c). 

Note that we have not displayed the contour plots of C(p) versus 6 /it and Jt for the 
Gaussian pulse magnetic fields, because they provide no extra information compared to 




' 5 10 15 20 



Figure 5. (Colour online) The contour plot of C{p) versus Jt and a/ ft with 9 = ~ and 
Jr = 0.1 (a), Jt = 0.2 (b) and Jr = 0.3 (c). Here the contour plots are for |11) initial 
state and include three positive Gaussian pulses centered at times T\ = 5, T 2 = 10 and 
T3 = 15, and we set J — 1 and (5 — 1. In these contour plots we have assumed twenty 
equidistant contours of concurrence between (black) and 1 (white). 



the ideal kick case; for the separable |11) initial state, C(p) is independent of 9 regardless 
of the value of the pulse width r. For the initial Bell state, the overall ^-dependence is 
very similar to that of the kick sequence control. 

We have also considered |00) separable, |$_) = ^(|11)-|00)) and |* ± ) = ^(|10)± 
1 1) ) Bell states as initial states and obtained some interesting dynamics. The dynamics 
of 1 00) initial state under control sequences considered in the present study are found 
to be exactly same as that of |11) initial state discussed above. From the propagators 
of Appendices and Eqs. (7|) and (8), one can easily show that |\I> + ) = J^(|10) + |01)) 
initial state is immune to the applied transverse field; its concurrence remains 1 for 
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(a) (b) 




\-\ ................ . . — ^ 

5 10 15 20 

m 

Figure 6. (Colour online) The contour plot of C(p) versus Jt and a/j3 with 6 = | 
and Jr = 0.1 (a), Jt = 0.2 (b) and Jt — 0.3 (c). Here the contour plots are for 
4g(| 11) + 1 00)) initial state and include three positive Gaussian pulses centered at 
times T\ — 5,T2 = 10 and T% = 15, and we set J = 1 and f3 = 1. In these contour 
plots we have assumed twenty equidistant contours of concurrence between (black) 
and 1 (white). 



every possible value of 9 and a/fi. Note that entanglement of this initial state can 
be manipulated by using a longitudinal control field as shown in Ref. [20]. Note that 
the interaction Hamiltonian, Hint ~ cos9(af + <rf) + sm.9(d~\ + erf) locally transforms 
Bell states among themselves for 9 = and 9 = tt/2. For example <r§. [$+) = 
and of. |$_) = |^_) as well as |$_) = |^+) and -io-f. [$+) = |^_). It is well 

known that entanglement is unaffected by local transformations, so one can easily show 
that the dynamics of |$+) at 9 = is the same as that of [$_) at 9 = ti/2 and 
concurrence remains 1 because the transverse fields transform the two initial states to 



Entanglement Dynamics of Kicked Qubits 



14 



\^+) = -^(1 10) + | OX) ) which is found to be unaffected by the control sequence as 
discussed above. Also |$ + ) at 9 = n/2 and |$_) at 9 = are transformed to \^-) 
and have similar entanglement dynamics; the concurrence for the initial state |^_) is 
independent of 9. At the intermediate values of 9, the local transformations mix different 
Bell states and it is not straightforward to unentangle them. All of them can be easily 
checked by using the propagators in the Appendices and the Eqs. ([TJ) and Q . The effect 
of transverse fast pulses on the entanglement dynamics of 1 10) or 1 1 ) separable states 
is qualitatively similar to the results for longitudinal fast pulse case [20], thus we have 
not covered them. 

It was pointed out that the observable non-local time-ordering effects in time are 
expected to be present if and only if the commutator of the total Hamiltonian Q at 



different times is nonzero, i.e, [H(t"), H{t')\ ^ 
show that the commutator [H(t"), H(t')} is equal to 



[H(t"),H(t')] = JD 



By using Eq. fllh, it is easy to 






a 


—a 





-a* 








a 


a* 








—a 





-a* 


a* 






xr 



Eq. (11) shows 



where D = ((B 2 (t') - B^if)) - (B 2 (f) - B 1 (t"))) and a = e~ ib 
that the commutator vanishes for the cases when J = and/or B\(t) = B2(t). By 



using the numerical solutions of Eq. (10) for the considered initial states and the 



propagators (A. 2), (B.2) and (C.2) in Eqs. ^ and (j8]), in the case of no qubit-qubit 
interaction (J = 0) as well as equal magnitude external fields (a = /3), it is easy to show 
that the concurrence for Bell state is equal to 1, while |11) state remains separable at all 
times and C{p) of both initial states is unaffected from the external kicks and Gaussian 
pulses. This situation can be also observed from the contour plots for a = j3 = 1 case. 
This finding indicates that ability to manipulate the entanglement in the system by 
using external fields is closely related to the nonlocal time-ordering effects. 



6. Conclusion 



We have investigated the possibility of creation and control of entanglement between 
two coupled qubits by using time-dependent external magnetic fields in x-y plane in 
the form of delta function kicks and Gaussian pulses for initially pure separable and 
maximally entangled states. Analytical (for kick sequence) and numerical (for Gaussian 
pulses) results presented and discussed in the paper indicate a number of interesting 
phenomena. 

Transverse fast pulses can be employed to create long lasting steady entanglement 
between two coupled qubits initially in a separable state |11) (or 1 00) ) . Furthermore, 
entanglement of such a state can be finely-tuned by changing the integrated magnetic 
strength of the external field at qubit positions, while the direction of transverse pulses 
are found not to effect the entanglement of this particular initial state. Similarly, 
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entanglement of a system initially in a Bell type state is found to be tuned with the same 
external control parameters; one can even destroy the initial entanglement and recreate 
this entanglement for this initial state with carefully designed pulse sequence and system 
parameters. We have also demonstrated the possibility of obtaining constant long-lived 
entanglement with desired magnitude by perturbing the qubits with fast pulse or pulses 
and by adjusting the system parameters to certain values for both initial states. 

The time-ordering effect defined by the commutator [H(ti), Hfo)] is found to be 
important in the manipulation and control of entanglement by fast pulses. It is found 
that in the case of no time ordering (i.e., [H(ti), Hfo)] = which holds for the cases 
J = and/or a = 0), the transverse fast pulses are found to be ineffective in creation 
and manipulation of entanglement between qubits initially in any type of pure states. 

We have also compared the effect of the external field being longitudinal or 
transverse and found that with a transverse control one can both create and manipulate 
entanglement, while the longitudinal fast pulses can only manipulate the entanglement 
of a state which is already nonzero [20] . Moreover, the transverse fast pulses enable one 
to manipulate, create and control entanglement for a number of initial states, while the 
manipulation of entanglement with the longitudinal fast pulses was found to be done 
only with a limited class of initial states [20] . 

Longitudinal and transverse control fields in the form of kick and Gaussian pulse 
sequences can be used in a proper sequence to create, control and destroy entanglement 
from arbitrary initial pure states of two coupled qubits. The same formalism can be 
applied to initially mixed states by using the time evolution matrix elements provided 
in Appendices to evolve the density matrix of the system as p(t) = U(t)p(0)W(t) which 
might be interesting in analyzing the time evolution of more-general-than entanglement 
type quantum correlations in mixed states, such as quantum discord [36] . 



Appendix A. Single positive kick 

Here, we consider two qubits whose states are strongly perturbed by external fields 
which may be expressed as a sudden kick at t — T. The time dependent magnetic 
fields on qubits 1 and 2 may be expressed as Bi(t) = aS(t — T) and B 2 {t) = (35(t — T), 
respectively, where a and /3 are called integrated magnetic strengths. For such a kick 
the integration over the time is trivial and the time evolution matrix in Eq. (|5| becomes 
as [201 E3]: 

U K (t) = e -Wo(t-T) e -if^:H int (t>)dt' e -iH T 7 (A1) 

with matrix elements 

f/ 11=e - lJt cos(|)cos =U U , 

1 <Jtf„4iJt ^ M 



U22 = 2 e~ lJt [e™ cos + cos ^-j j = U 33 , 
U* = - le~ ut ( e iUt cos ( ^] - cos f ? I ) = U 



2-4 ~2 \ \ •_> ./ ~ V 2 ' ' 
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U 12 = \ % e-^e~* (e 4 ^sin (|) - sin (|)) = e" 2 ^, 
E/ 13 = - -ze^e"* (e 4iJT sin ( A )+sin (~\) = e~ 2ie U, 



i q — te e e sin — ~r sin — — e (-742, 

2 V V2/ \2// 

iJt-2iO- l' a \ • ( P \ -4i0; 



[/14 = - e^e-^sin sin I ^ J = e^U^ 
U 2l = - C <J <«^e* (cos sin (|) sin(£) + i cos (|) sin (A cos(o) 

f/ 24 = - (e 2 ^ sin (|) + sin (^)) = e~™U 3U (A.2) 

where A = (a — (3),Q = (a + j3) and £ = 2J(t — T). It should be noted that the 



propagator given by Eq. (A.l) is valid only at times t > T. 



Appendix B. Two positive kicks 

The next example is the positive-positive kick sequence applied at times t = T\ and 
t = T 2 , namely B^t) = a(5(t - T t ) + 5(t - T 2 )) and B 2 (t) = (5{5{t - T x ) + 5(t - T 2 )). 



Following the procedure given in Eq. (A.l), one obtains the time evolution matrix at 
times t > T 2 as [SB [23]: 

JJ K (t) = Q-iHoit-Ti) e -i Ull-l H mt {t')dt' e _^ o(T2 _ Tl ) e -i H int (t')dt' e _iH Tx ^ 

(B.l) 

with matrix elements 

U u = \e~ iJt (l + e iiJT (cos(A) - 1) + cos(A) + 2 cos(fi)) = U u , 

U 22 = ^e~ iJt (e 4iJt + e 4iJ( *" T) (cos(A) - 1) + e 4lJt cos(A) + 2cos(fi)) = U 33 , 

U 23 = - ^e~ iJt (e 4iJ * + e 4iJ{t ~ T) (cos(A) - 1) + e 4iJt cos(A) - 2 cos(fi)) = U 32 , 

U 12 = Ue- iJt e~ ie (e 6iJT cos(2 JT) sin(A) - sin(fi)) = e~ 2i9 U A3 , 

U 13 = - l -ie- lJt e~ ie (e 6i JT cos(2 JT) sin(A) + sin(fi)) = e" 2 ^, 

[/ 14 = \e- iJt e-™ (e 4UT - 1 - (1 + e 4iJT ) cos(A) + 2cos(fi)) = e~ 4ie U 41 , 

U 2 i = \%e~ iH e ie ( e 4iJ(i " 2T) (l + e 4iJT ) sin(A) - 2sin(fi)) = e 2ie f/ 34 , 

U 24 = - ^e- m e- w (e 4tJ (*- 2T )(l + e 4lJT ) sin(A) + 2 sin(fi)) = e^E/si, (B.2) 

where A = (a — 0) and Q = (a + f3). Here, we have assumed equally distanced kicks 
applied at times T\ = T and T 2 = 2T. 
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The final example is the sequence of three positive kicks applied at times t = T, t = T 2 , 

3 3 

and t = T3 namely, B\(t) = a8(t — Ti) and B^it) = ^ 08(t — Ti). Following the 

i=l i=l 

procedure given in Eq. ( |A.l ), one can obtain the time evolution matrix at times t > T 3 
as EDI: 



with matrix elements for T\ = T, T2 = IT and T3 = 3T, 

U u = ~e~ iJt ((3 - 2e 4UT - e 8UT ) cos (~) + (1 + e 4UT ) 2 cos + 4 cos ) 



U, 



44- 



E/22 = \e- m (cos + e 4 ^'"^ cos (~) ((1 + cos(4JT)) cos(A) + i sin(4JT) - 1] 



U33 



U 23 = \e~ m (cos - e 4iJ( *" T) cos ((1 + cos(4JT)) cos(A) + i sin(4JT) - 1) 

= U 32 , 



C/ 12 = -i e - JJt e-^ ^e SUJ (cos(4JT) + 2 cos(2 JT) 2 cos(A)) sin (-) - sin _ ? y 



43) 

1 . _»• 7+ _»o / a.' TT / /. Tm \ „ /„ Tm \ o / » \ \ . /A\ . / 3f2 



[/ 13 = - - ze - iJt e- ie (e 8iJT (cos(4 JT) + 2 cos(2 JT) 2 cos(A)) sin + sin 



2 



1 - iJt e- 2iB ( (2e 4iJT + e 8lJT - 3) cos ( - (1 + e 4iJT ) 2 cos (^) + 4 cos 



M2 

£/ 14 = ^ e - lJ V M " ( (2r " Jl +< > '-' 1 - 3)cc.s( = ] - 1] + f )-c-(,s I — j - I cos 
— e C41, 

l/ 21 = Ue- m e lf> ( e 4iJ (*- 2T ) (cos(4 JT) + 2 cos(2 JT) 2 cos(A)) sin f^) - sin 



e 2W U- 



34- 



f/ 24 = - ^e- iJt e- ie fe 4iJ ( i - 2T Hcos(4JT) + 2 cos(2 JT) 2 cos(A)) sin (^) + sin ^""^ 



2 J, 

= e~ 2l6 U 31 , (C.2) 
where A = (a — 0) and Q — (a + 0). 
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